ANGULAR DEPENDENT VORTEX DYNAMICS 

IN SUPERCONDUCTORS WITH COLUMNAR DEFECTS 



o 

^ ■ LEONARDO CIVALE 

Superconductivity Technology Center, 
^ I Los Alamos National Laboratory, Los Alamos, NM 87545, USA. 

< 

^ ; ALEJANDRO V. SILHANEK 

CN ! Laboratorium voor Vaste-Stoffysica en Magnetisme, 

K. U. Leuven, Celestijnenlaan 200D, B-3001 Leuven, Belgium. 

g : GABRIELA PASQUINI 

^ ' Laboratorio de Bajas Temperaturas, Departamento de Fisica, 

Q^' Facultad de Ciencias Exactas y Naturales, 

^ , Universidad de Buenos Aires, Buenos Aires, Argentina. 

a 

o 

CJ : I. INTRODUCTION 

Motivation 

> 

in 

O ! The introduction of aligned columnar defects (CD) in high temperature superconductors 

(HTSC) was initially motivated|l|] by the technological objective of increasing the discour- 
^ ■ agingly low vortex pinning of the early samples of these compounds. Initial results clearly 
! showed^, 0] that CD were extremely effective in improving the critical current density (Jc) 
(particularly at high temperatures and magnetic fields) and in enhancing the irreversibility 
c5 ■ field (Hirr). In addition, it soon became clear [3, 4] that good quality single crystals with CD 
were excellent "model systems" to explore the complexities of vortex dynamics in HTSC. 
The combination of both motivations produced a very large interest on the topic|5:], which 
Pi I remains active after more than a decade of research, 
g ' Vortex pinning in type II superconductors arises 0, 1^ from the presence of crystalline 

defects, which produce a spatially inhomogeneous superconducting order parameter '?/'(r), 
thus inducing a position dependent vortex energy. Examples of those defects are dislocations, 
^ I chemical doping, vacancies, non superconducting particles, voids, etc. The calculation of 
d • the expected Jc as a function of temperature (T) and magnetic induction (B) is a very 
complex problem^, |^, as it depends nontrivially on a large number of variables such as 
the nature of the defects, their size, shape and density, as well as on the temperature and 
magnetic field dependent elastic properties of the vortex array, which in turn are determined 
by the superconducting parameters of the material. If the material and/or the defects are 
anisotropic, pinning will additionally depend on the orientation of the applied magnetic field 
H. 

A particularly important type of anisotropic pinning is that produced by extended parallel 
pinning potentials, commonly known as correlated pinning^. In HTSC, common examples 
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of correlated pinning structures are the CD, the twin boundaries [tI. IM IqI ll^ 11 1, and the lay- 
ered crystallographic structure of these materials that produces the intrinsic pinningjl^.[l3l|. 
Although each one has its distinctive characteristics (CD are linear while the other two are 
planar, intrinsic pinning potentials are spatially periodic while the separation between CD 
and twins is random), they share several important similarities in their pinning properties, 
particularly in regards to the vortex structures that occur0, |^ as a function of the orienta- 
tion of H. In all cases, pinning is maximum when the vortices are parallel to the defects. 
They all exhibit a lock-in phase: when the angle between H and the defects is smaller 
than a certain lock-in angle ifL-, the vortex cores remain locked into the pinning potentials. 
For if > ifL, vortices form staircases, with core segments locked into the defects connected 
by weakly pinned kinks. 

Different mechanisms can produce aligned linear defects capable of pinning vortices in 
a superconductor. In YBa2Cu307 (YBCO), for instance, pinning by dislocations 3 [3| 
and by the lines occurring at the intersection of two perpendicular families of twins la] has 
been explored. However, we will restrict the term "columnar defect" to those produced by 
irradiation with very energetic heavy- ions. Several features set them apart from the rest. 
First, they can produce a much larger Jc than any other types of linear defects (as discussed 
in section I.C. below, this is mainly due to the fact that CD have the appropriate diameter). 
Second, their density can be easily controlled. Third, they can be introduced at an arbitrary 
crystalline orientation. This is a unique advantage of the CD as a tool to explore vortex 
dynamics, as it allows to deconvolute the features arising from their uniaxial pinning from 
those due to anisotropy and other correlated pinning mechanisms. We also note that by 
"aligned" CD we do not mean that all the tracks are perfectly parallel, as will be discussed 
in the next three subsections. We just use this term as opposite to CD pur pos ely generated 
with high angular dispersion 17|, |l8| (splay) or even randomly oriented [l9|, |20|. 

The purpose of this review is to discuss the angular dependence of the vortex dynamics in 
type II superconductors with aligned CD. The problem is rather complex, and some aspects 
still require further investigation. Our main goal is to identify and characterize the various 
angular regimes, both statically and dynamically. To that end we will use and combine 
results obtained by dc magnetization and ac susceptibility measurements. We will show that 
both the crystalline anisotropy and geometrical anisotropy (sample shape) play an important 
role. Most of our studies have been performed on high quality twinned YBCO single crystals 



prepared using a flux growth technique j21| . In this compound, the combined effects of the 
three types of correlated pinning mechanisms and the anisotropy produce a rich variety of 
vortex structures and dynamic responses. Complementary studies were performed on high 
quahty single crystals 22] of the less anisotropic conventional superconductor NbSe2, with 
Tc = 7.3 K, coherence length ^ = 77 A, and penetration depth A = 690 A, and having the c 
axis parallel to the thinnest dimension. Magnetization and susceptibility measurements were 
performed at the Centro Atomico Bariloche and Centro Atomico Constituyentes respectively, 
both dependent on the Atomic Energy Commission of Argentina. 

This review is organized as follows: In the remaining of section I, we briefly describe 
the generation of CD by heavy-ion irradiation (section LB), and we summarize the basic 
vortex dynamics for superconductors with aligned CD (the case H ]] CD is discussed in 
I.e., and the angular dependence in I.D.). In section II we present an overview of the 
temperature, field and angular dependence of the irreversible magnetization in samples with 
columnar defects and analyze the influence of correlated defects on the pinning properties of 
vortex lines. First (in section II. A.) we discuss the origin of the irreversible magnetization 
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in superconductors and introduce our measurement procedures. In section II. B. we focus on 
the angular regimes associated with different vortex configurations as the external field is 
tilted away from the CD orientation. In section II. C. we show that the sample geometry and 
the material anisotropy are competing mechanisms that determine the vortex orientation at 
low fields. In section III we use our ac susceptibility measurements to explore the vortex 
dynamics in YBCO with CD in the vicinity of the solid-liquid boundary. In section III. A. 
we introduce the experimental procedure. In section III.B. we study the angular dependence 
of the ac response with the scope to identify the various angular regimes and to investigate 
which is the main source of pinning at different orientations of a low DC field. In section 
III.C. we explore the dynamic response of vortices aligned with the CD over a wide range of 
current densities and excitation frequencies, and we construct the dynamic phase diagram 
of the system in the temperature vs ac field plane. In section III.D. we extend that study 
to different angles between the applied field and the tracks. 

B. The generation of CD by heavy-ion irradiation 

The generation of CD by irradiation with very energetic heavy ions has been extensively 
discussed in the literature (for a review see for example ref. and references therein). Here 
we will only present a very brief summary of the basic concepts. 

Particle irradiation is a standard way to introduce defects in a solid. Incident particles 
transfer energy to a solid psi. by direct collisions with lattice atoms (nuclear or 

nonionizing energy loss, Sn), and by ionization or electronic excitations (electronic or ionizing 
energy loss, Se)- In the case of irradiation with electrons, protons and light ions with energies 
up to a few MeV, almost all the energy is transferred via nonionizing energy loss [2^. This 
results in defect formation through displacement of the primary recoil atom, and through the 
cascade of collisions that this atom produces if its energy is high enough |25|. The size of the 
damaged region increases with the primary recoil energy transfer. Electron irradiation (up 
to a few MeV) can only produce point defects (vacancy- interstitial pairs or Frenkel pairs), 
while protons and light ions generate defects up to several tens of A in size. Fast neutrons 
also generate cascade defects, with diameters in the range of 50 to 100 A. In all cases, these 
defects are randomly distributed. 

As heavier and more energetic ions are considered, Se grows and the electronic energy 
loss becomes dominant [24] . For instance, for our typical irradiations using ~ 300 MeV 
Au, Se ~ loos'™, while for ion energies in the GeV range the relation can be as large j26l| 
as Se ~ 20005'„. In this limit each incident ion forms a cylinder of amorphyzed material 
along its path, around 40 to 100 A in diameter, known as a latent track. The heavier 
and more energetic the ion, the more continuous, homogeneous and parallel are the tracks. 
In fact, as Se increases the defects evolve jl^, ll^l from isolated uncorrelated spheres to 
aligned spheres, then to elongated aligned but disconnected defects, which then coalesce 
into discontinuous inhomogeneous tracks, which eventually become continuous and finally 
homogeneous. However, as will be pointed out below, perfectly aligned and homogeneous 
defects are not necessarily the best suited for high vortex pinning. 

The Se values separating the various regimes depend on the material. The thresholds are 
lowest in insulators, while in the other extreme it is impossible to create continuous tracks 
in good metals. For the case of YBCO, irradiation with 173 MeV Xe {Se ~ 1.2 KeV/A), for 
instance, does not generate defects capable to produce uniaxial pinning. In this material, 
the threshold for the formation of continuous inhomogeneous tracks isQ ~ 2 KeV/A, 
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somewhat below the 2.7 KeV/A that occurs at the entry surface for the 580 MeV Sn"^^^ 
used[l| in the initial studies. Heavier ions such as Au and Pb in the GeV energy range, 
which have Se as highjis^ as 3.5 to 4.5 KeV/A, result in homogeneous tracks. All these Se 
values are extremely large, so incident ions lose energy very fast as they penetrate deeper 
into the target. As they do, at some point Se starts to decrease, so the morphology of 
every track changes as a function of depth, until the ion is eventually stopped. Typical 
range of continuous tracks formation is about 10 fim for our 300 MeV Au, and can be as 
large as ~ 100 /im for energies in the several GeV range. As ions penetrate deeper into 
the material, the angular dispersion of the tracks also increases ji^. This is mainly due 
to Rutherford scattering of the incident ions, associated to the small but nonzero nuclear 
energy loss. A completely different approach based on fission tracks can be used[l^, li^l to 
generate randomly oriented CD. 

Irradiations reported in this review were performed at room temperature at the TANDAR 
accelerator facility. A group of samples was irradiated with 291 MeV An'^'^+ ions and another 
one with 220 MeV Sn^^+ ions. In both cases, the energy deposition rate is greater than the 
threshold for continuous latent track formation in the first 9 — 10 fim. So we only used 
crystals of thickness < 10 fim. 

C. Vortex dynamics in superconductors with CD: basic results for H || CD 

In this subsection we present a brief summary of the basic results of vortex dynamics 
when the applied magnetic field is parallel to the CD. More detailed theoretical analyses 
can be found in refs. 3 and 4. 

The basic reason why CD are very effective pinning centers when the vortex direction 
is parallel to them is quite simple: each CD can in principle confine the whole length of 
the vortex core, thus producing the maximum possible pinning energy without any cost in 
elastic energy due to vortex bending. This contrasts with the case of randomly distributed 
defects, where only a fraction of the core length is pinned, and additionally vortices must 
zig-zag among defects in order to get pinned, thus paying an additional price in terms of 
vortex elastic energy. 

Let's start with a rough estimate. At low temperatures and for an isolated vortex {i.e., 
at low enough fields) the pinning energy per unit length of a CD (assuming ip = inside it) 
is Ep ~ (iJ^/87r)7r6^, where He is the thermodynamic critical field and b is the smallest of the 
defect radius and the core radius \/2^ (here ^ is the superconducting coherence length). 
The pinning force per unit length will be of the order of Fp ~ Ep/^, and the critical current 
density Jc = cFp/^Q, where $o is the flux quantum and c is the speed of light. Clearly, 
the most convenient case is ro ^ v^^, so the highest Fp is achieved without removing 
more superconducting material than necessary. In this case we obtain Jc ~ </o/2, where 
Jo = cHc/S^/GttX is the depairing current and A is the superconducting penetration depth 
for H II c axis. A more exact calculation gives0, 0] Jc ~ (3-\/3/4\/2) Jq ~ Jq, showing 
that, at least in principle, CD can produce a Jc pretty close to the maximum theoretically 
achievable value. It is a fortunate fact that the CD produced by heavy-ion irradiation have a 
radius ~ 20 — 50 A that matches very well with the vortex core radius ~ 25 A of HTSC. 
On the other hand, the reality is that the largest Jc's measured^, at low temperatures in 
YBCO crystals with CD are ~ 4 — 5 x 10^ A/cm^. Although these values are very high, 
they are far from the ideal limit ~ Jo{T = 0) ~ 3 x 10^ A/cm^. This may be due to a 
variety of factors related to the morphology of the tracks, such as its inhomogeneities, the 
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degree of depression of ib, the sharpness of the normal/superconducting boundary, etc. It is 
thus useful to introduce j3(J| a dimensionless efficiency factor 77 < 1 in the pinning estimates, 
which accounts for all the uncertainties mentioned. 

Several additional factors reduce Sp (and thus Jc) as temperature increases 0, 0|- First, 
HciT) decreases with T. Second, as ^(T) increases with T a crossover from r^) > \/2^ to 
rn < -\/2'C may occurs. The third and most important factor is due to thermal fluctuations 
and is known as entropic smearing. The transverse localization length of the vortex cores, 
which is determined by the root mean square amplitude of the thermal fluctuations, grows 
as ksT, thus reducing the vortex free energy F = U — TS due to the increase of the entropy 
S arising from the delocalization. As a consequence, the effective potential well associated to 
the CD becomes "shallow" and the pinning energy decreases by a factor f{x), where f{x) ~ 
exp(— x) for X > 1, according to the long range nature of the CD-vortex interactionj^. 
Here x = TjTdp and Tdp is called the depinning temperature. In the very useful analogy to 
2D-bosonsp|, this corresponds to the decrease in binding energy of a quantum well due to 
the zero-point fluctuations. The entropic smearing effect has been experimentally observed 
in YBCO, and it was found (sol Isi ] that Tdp ~ 40 K. The combination of all the factors 
mentioned above produce an effective pinning energy per unit length 

^p(T)=r/|ln(l + ^) x/(a:) (1) 

where So{T) = ($o/47rA)^ is the vortex energy scale, and the logarithmic factor provides a 
convenient interpolation P| between the td > and rjj < cases. 

As the vortex density increases, more CD will be occupied. The natural characteristic 
field in the problem is[l| the matching field = n^Q, where n is the areal density of 
CD. For B ^ B^ there are many more vortices than CD, so all the CD are occupied and 
only a small fraction ~ B/B<j, of the flux lines are actually trapped, the rest of them being 
pinned only as a result of vortex-vortex interactions. We thus expect Jc to decrease fast 
with B in this regime, as indeed observed 0, |^, %0]. When B <^ Bq, each flux line will 
be confined into a CD. However, as B increases it may occur that, although there is a 
CD available, it will be energetically unfavorable for the vortex to occupy it because it is 
too close to another vortex, so the gain in pinning energy will be overcompensated by the 
increase in the elastic energy associated to the vortex-vortex repulsion. The characteristic 
field where this effect sets in is called the accommodation field, Ba{T) oc {ep{T)/eo) Bi},, 
which represents the crossover between single vortex pinning where vortex interactions are 
negligible and collective pinning where they become significant 0, 0, |3^ . At low temperatures 
Ba approaches B^, but above Tdp it decreases due to the reduction in Sp(T). Above a 
temperature Tdi > Tdp the localization length becomes larger than the average distance 
between tracks, so vortices become collectively pinned by several tracks and Ba{T) decreases 
faster. We have estimatedjlj that in YBCO Tdi ~ 70K. 

Thermal fluctuations are also responsible for the large flux creep in HTSC. In the presence 
of CD and below Ba{T), initial relaxation proceeds via nucleation and expansion of half 
loopsjl, 0, H^. This is a glassy regime, in the sense that the current density dependent 
activation energy U{J) increases as with decreasing current density. Below a certain 
J, a crossover to a creep regime dominated by double kinks occurs j33|. These excitations 
involve the transfer of a flux line from one CD to another one, and once a double kink 
is formed its further expansion is energetically favorable, thus relaxation in this regime is 
fast and nonglassy. Below a still lower J, the expansion of double kinks is precluded by 
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the dispersion in the pinning energy of the CD. Relaxation now evolves through another 
type of excitations called superkinks. This regime, analogous to the variable range hopping 
conductivity in doped semiconductors, is again glassy, with U{J) diverging as J~^/^. For 
B > Ba{T) intervortex interactions are relevant at all current densities, and the relaxation 
takes place via collective excitations involving vortex bundles. A description of the relaxation 
regimes in that case can be found in Ref. 0- 

An interesting implication of the above analysis is that identical and perfectly parallel 
CD are not the ideal configuration to sustain large persistent current densities, because 



. In fact, it has been 
that larger persistent 



they allow for the fast, nonglassy relaxation of double kinks y. 
theoretically predicted jSJ], and experimentally demonstrated [l 
currents can be obtained in the case of CD having a small angular dispersion or splay. The 
case of very large angular splay is also interesting and has been explored by generating 
randomly oriented fission tracks |l9l. 2ol 



A related phenomenon is the matching effect. When B = B^, the available phase space 
for transfer of vortex cores between two CD should be strongly limited. Thus, creep rate 
should be reduced and the field dependence of the persistent current should exhibit a peak 
at -B$. However, this matching effect is generally not observed^, |^. The reason is that 
usually the energy and angle dispersion of the CD reduces the creep at all fields 3^ , thus 
masking this phenomenon. Matching effects have only been observed [ssl. l36l | in very thin 
samples containing CD with very low energy and angle dispersion. 



D. Angular dependence: Bose glass picture 

Our studies of the angular dependence of vortex pinning in HTSC with CD have revealed 
a richer variety of phenomena and pinning regimes than originally expected. As we will 
extensively discuss in the following sections, this is due to a number of additional ingredients 
such as influence of other types of correlated disorder, energy dispersion of the CD and 
misalignments between B and H due to both mass anisotropy and sample shape. But 
first, in this subsection we briefly summarize the expectations within the basic Bose Glass 
scenarioj^, 0| for a single family of identical and perfectly parallel CD. In principle, this 
picture also applies (with minor differences) to twins and intrinsic pinning. 

First, when the angle between the applied field H and the correlated defects is smaller 
than a certain ipi, it is energetically favorable for the vortices to remain locked into the 
defects. This is related to the concept of transverse Meissner effect. The lock-in angle is 

where Si is the vortex line tension. In isotropic superconductors Si = Eq In k, but in the 
anisotropic case the situation is rather complex. In our experiments the appropriate line 
tension is that corresponding to in-plane deformations (see pages 1163-1164 in Ref. 0), 
El = (£^£o/^(0)) In K, where the mass anisotropy e = niab/mc -C 1 and e^(0) = cos^(9) + 
e^sin^(G). Note that both ei and Ep decrease with T. 

For tilt angles larger than (pi and smaller than a trapping angle, (pr = ^/'^^p/^h vortices 
form staircases with segments pinned into different defects and connected by unpinned or 
weakly pinned kinks. Beyond y^r, vortices become straight lines aligned with the applied 
field and unaffected by the correlated nature of the pinning potential. 
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Experimentally, the determination of bulk vortex structures for applied fields tilted with 
respect to the correlated pinning potential is very difficult. Most of the existing imaging 
methods that permit direct observation of the vortices can only detect the flux lines at 
the sample surface js^]- Recent studies have shown that it is possible to resolve individual 
vortices inside a superconductor with columnar defects by means of Lorentz microscopy 
or interference microscopy 3^ . Unfortunately, however, the low penetration power of the 
electron beam imposes a strong restriction in the maximum thickness of the sample of about 
0.5 /im. Thus, in general, in order to obtain information from the 3D nature of the vortex 
structure we still have to rely on indirect techniques like dc-magnetization, ac-susceptibility 
or transport measurements [ll.l7ll8l.l9l ll^ 11, 4^, 41, 4^, 4^, i^]. In this context, it is useful to 
introduce CD tilted off the crystallographic c axis, as this helps to discriminate the pinning 
produced by them from that due to twin boundaries and from anisotropy effects. 



II. DC MAGNETIZATION: ANGULAR DEPENDENT VORTEX 
CONFIGURATIONS IN THE SOLID PHASE 

A. Introduction to dc-Magnetization Measurements 

The widely used dc magnetization measurements are a powerful tool to explore the ir- 
reversible response of the vortex system deep into the vortex-solid phases, where electro- 
transport measurements and ac-susceptibility hardly access. If a superconductor is cooled 
down in a zero applied field (ZFC experiment) and then an external field H > Hd is ap- 
plied, fiux-bearing vortices enter through the sample's border until their motion is arrested 
by pinning centers. As a consequence, the system achieves an inhomogeneous flux distribu- 
tion with a higher density of vortices near the border that progressively decreases toward 
the center of the sample. The spatial variation of the locally averaged field B(r) gives rise 
to supercurrents J in the sample that, in the stationary state, accommodate to be exactly 
the critical current Jc everywhere. This is the so-called Critical State regi nie p . lisj] . 

The simplest version of the critical state is the isotropic Bean model j45|, where Jc is 
assumed to be field independent and isotropic, so the current density throughout the sample 
(in the fully penetrated state) is necessarily uniform, and the associated irreversible magne- 
tization Mi is proportional to Jc via a geometrical factor. For instance, for a disk of radius R 
with H parallel to its normal n we have Mi = {JcR/3c) n. However, even in this simple case 
the calculation of the field profile is rather complex, and has been the subject of extensive 



modelling and numerical analvsis|46l. 1471 . |48[. The situation is further complicated when 
realistic cases of field dependent and anisotropic critical currents are considered |49l. Isoj] . 

If H is tilted with respect to n, then Mi has two components. My and M^, parallel 
and perpendicular to n respectively (see Fig. that in general depend on the various 
components of Jc and on the sample geometry. In very thin samples (which is our case), 
since the non-equilibrium screening currents are constrained to flow parallel to the sample 
surface, there is a large angular range of applied fields 51, 5^ in which Mi points almost 



parallel to n. For an infinite strip of an isotropic superconductor, Zhukov et a/.^] showed 
that Mi remains almost locked to n as long as the angle Qh between H and n remains 
smaller than Be = arctant/"^. Here u = 6/w, where 6 and w are the thickness and width of 
the strip respectively, with 5 -C w. For this particular geometry. 
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Mx_ = — JctanQH Mn = | 1 
be 4c 

and the angle 6m between Mi and n is given by, 



tan^e^ if eH<Qc 



tan 



M 



M_L 



2z/^ tan Qh 
2 tan2 e 



H 



-v^ tan Qii 



(3) 



(4) 



Typically we deal with samples of dimensions 5 ~ 15 /im and w ~ 400 /im, so z/ ~ 0.04 
and Gc = 87°, i.e. M\\ ^ M_l and the geometrical factor changes only 0.1% for < Be- 
On the other hand, for B/^ > Be the irreversible moment vector suddenly rotates towards 
the sample plane. 

It must be noted that the total field in the sample results from the applied field plus the 
self field Hgf generated by the screening currents. At the center of a thin disk H^f ~ J^t 
(rather than ~ X.i?)and the shielding currents also create radial fields of ~ Hsf/2 on the 
disk surface [iil. WA lisjl . For H < H^f self field effects dominate, thus Jc{H) cannot be easily 
extracted from Mi{H). We constrain our analysis to the regime H ^ H^f in order to obtain 
reliable estimations using the Bean model. For instance, at T = 60 K for YBCO crystals of 
6 ~ 15/im with CD, ~ 10^ A/cm^ and hence Hsf ~ 150 G. 

In spite of the shortcomings of the Bean model and the necessity of extend this description 
to account for anisotropics and field dependencies of Jc, several experiments have shown that 
the behavior of the irreversible magnetization follows accurately this model 5^ 5^ justifying 
its use as a reliable way to determine the critical current of the system 5^. 



The magnetization measurements presented in this section were conducted on a SQUID- 
based magnetometer Quantum Design MPMS-5S equipped with two sets of detectors, which 
allows us to record both the longitudinal Ml{H) and the transverse Mt{H) components 
(parallel and perpendicular to H, respectively) of the total magnetization vector M. The 
crystals can be rotated in situ around an axis perpendicular to H (see Fig. P), and they are 
carefully aligned with the rotating axis normal to the irradiation plane, in such a way that 
the condition H 11 tracks could be achieved within ~ 1°. 




FIG. 1: Sketch showing the orientation of the irreversible magnetization Mj and applied field H 
with respect to the sample normal n. 
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Our usual measurement procedure [56|, |57| consists in a ZFC of the sample from above 



Tc down to the work temperature T at a given Qh- Then, both components Ml{H) and 
Mt{H) of a magnetization loop are recorded (at fixed T and Qh)- After that, the sample is 
rotated to a new desired angle, warmed up above and cooled down in zero field to start 
a new run. The total magnetization vector so obtained is M = Mi + Mgg, where Mgg is 
the equilibrium (reversible) magnetization in the mixed state of the superconductor, plus 
any contribution from the sample holder. We use the widths of the hyterisis AMl{H) and 
AMriH) to calculate the modulus Mi = /\Ml{H) + AM^{H) and the angle Qm = 
Atan{AMT{H) / AMl{H)) of Mi(H). Note that in this way M^q automatically cancels out 
of the calculation, so it is not necessary to determine it. 

We have also developed an alternative method which allows us to obtain directly 
MiiQn) by rotating the sample at fixed H and T. This procedure is much faster, and its 
principal advantage is that a finer grid can be easily obtained in the angular range of inter- 
est, thus permitting the exploration of various regimes with significantly improved angular 
resolution. The complication is that the M^q is not automatically cancelled, thus it must 
be subtracted (usually, however, in HTSC the strong pinning determines that M, ^ Meq). 
We have solved that complication by appropriate data analysis^S^, and we have been able 
to obtain very good agreement between data from sample rotations and from mag netization 
loops. A rotation at fixed H is to some extent analogous to a hysteresis loopj59|. Rotating 
the sample normal, fi, towards H increases if^, which is roughly equivalent to increasing 
H at Qh = 0°, moving along the lower (diamagnetic) branch of the loop. Decreasing H± 
(either by rotating n away from H or by crossing the H || c condition), is equivalent to 
reversing the field sweep, thus producing a switch to the other branch of the loop. This is 
a useful analogy for the analysis of the rotations, although it should not be pushed too far. 

Finally, it is important to mention that, due to the large influence of thermal fluctuations 
on the vortex dynamics in these HTSC compounds, the persistent current J determined 
through dc magnetization measurements in the typical time scale of SQUID magnetometers 
(~ 30 sec.) is much smaller than the "true" critical current density Jc- 



B. Angular regimes 



1. Lock-in Phase 

Figure |21 shows typical Mi versus Qh curves at T = 60 K for a YBCO single crystal with 
CD at an angle Qd = 32° away from the c axisj^. As discussed above, the geometrical 
factor between Mj and J is almost constant for < 0c ~ 87° in our crystal, thus the 
vertical axis in Fig. |21is directly proportional to the persistent J over almost all the angular 
range. The uniaxial nature of the pinning potential is clearly manifested as an asymmetric 
angular response M((Bh) ^ Mi{—QH)- At high fields {H > 10 kOe) we observe a large 
peak in the direction of the tracks Qmax — Qd = 32°. At lower fields the peak becomes 
broader and transforms into a plateau (the angular range where Mj(Gf/) ~ const.) as well 
as it progressively shifts away from the tracks in the direction of the c axis {Qmax < ©d)- 
Hereafter in this section we will focus in the study of the plateau and defer the discussion 
of the shift to the next section. 

The inset of Figure |21 shows a zoom in of the data of the main panel for H = 4 kOe, 
where the plateau is clearly observed. The constancy of Mj (Qh) indicates that the pinning 
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0„ (deg) 



FIG. 2: Irreversible magnetization Mi[H) as a function of Qh for several fields, at T = 60 K. Inset: 
zoom in of the data of the main panel for H = 4 kOe around the plateau region. 

energy remains constant and equal to the value at the alignment condition = Qd- 
This behavior is a fingerprint of the lock-in phasejs^. Moreover, the decrease of Mj at the 
edges of the plateau is quite sharp, a result consistent with the appearance of kinks, which 
not only reduce Jc but also produce a faster relaxation. A few examples of the observed 
plateaus for a ErBa2Cu307 (ErBCO) single crystal are shown in Figure El for several fields 
and temperatures (some curves have been translated vertically for clarity). 

The plateau represents the angular range of applied field over which it is energetically 
convenient for the vortices to remain locked into the columnar defects, therefore its angular 
width is twice the lock-in angle (pi- The results for all measurable plateau widths are plotted 
as a function of in Figure 0] This figure clearly demonstrates the linear dependence 
of on over the whole range of temperature and field of our study (as evidenced by 
the solid lines which are the best linear fits to the data), in agreement with the Bose-glass 
prediction [Eq. (j2))]. 

To analyze whether Eq. Q provides a satisfactory description of the temperature de- 
pendence of the lock-in effect, the slopes a(T) = dipi/ d{H~^) of the linear fits showed in 
Fig. lUare plotted in Fig. El (solid symbols). As expected, a(T) decreases with increasing T, 
reflecting the fact that the lock-in angle at flxed H decreases with T due to the reduction 
of both, the line tension and the pinning energy. From eqs. (Q) and (j21) evaluated at the 
track's direction = O/? we obtain 

where f{x) ~ exp(— x) for a; > 1. 

We can now flt the experimentally determined ot{T) using Eq. To that end we flx 
the reasonably well known superconducting parameters of the material e ~ 1/5; Iuk 4 



D 



(5) 
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and ^ = 15A j^X — t (where t = T/Tc). We also assume the usual two- fluid temperature 
dependence A(T) = Aj;,/2a/1 — t^, where Xl is the zero-temperature London penetration 
depth. The free parameters are then T^p and the combination Xl/v^^^- The best fit, shown 
in Fig. as a sohd line, yields Xl/i]^^^ = 360 A and T^p = 30 K. 

The obtained T^p is smaller, but still reasonably similar to the value ~ 40 K reported 



for several YBCO crystals using a completely different experimental method (30|, |31|. This 
low Tdp (well below the initial theoretical expectations) indicates that the efficiency factor 
1] is rather small, what is also consistent with the less than optimum Jc observed here and 
in several previous studies. For low matching fields as that used in the present work, it was 
estimated^ that r/ ~ 0.2 - 0.25. 

The exact value of r] has little influence in our estimate of Xl, as it only appears as t]^^^. 
For 7] = 0.2 and rj = 1 we get Xl = 250 A and 360 A respectively, a factor of 4 to 5 smaller 



than the accepted value A/, ~ 1400 A. Zhukov et a/.[ll| had reported a similar discrepancy 
when studying the lock-in effect by both CD and twin boundaries in YBCO. Similarly, as we 
will show in the next section, we also find that the misalignment between B and H at low 
fields can be described satisfactorily using a A^ smaller than the accepted value^]. Thus, 
this numerical discrepancy appears to be a common result associated to the study of angular 
dependencies in YBCO-type superconductors with correlated disorder that deserves further 
analysis. 

Finally, it is relevant to note that Eq. (0) was derived for the single vortex pinning 
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FIG. 4: Plateau widths ipL vs H for several temperatures. The straight lines are fits according 
to Eq. @. 
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FIG. 5: Slopes Q(r) = d^L/d{H-^) of the linear fits of Fig. H 



regime, which occurs below the temperature dependent accommodation field Ba{T) < Bq, 
(section I.C), while a large fraction of the data shown in Fig. |3]lies above this line, in the 
collective pinning regime. Unfortunately, to our knowledge there is no available expression 
for ifL{H,T) in the collective regime. Blatter et a/.|4:J only argued that collective effects 
should result in a reduction of the lock-in angle. The experimental fact is that Eq. Q 
satisfactorily describes both the temperature and field dependence of ^9^. This suggests 
that, at least to a first approximation, collective effects in the range of our measurements 
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simply result in a different prefactor in Eq. 
2. Trapping Angle 

As was mentioned in Section I.D., it is expected that for \Qb — ©dI > fL vortices form 
staircases. The question that arises here is whether we are able to observe any clue of 
the trapping angle ifT from the angular dependence of the irreversible magnetization. For 
clarity lets concentrate the analysis on high enough fields such that no plateaus or peak 
shifts appear. 

First, we have to notice that for > ©max, the asymmetry M, {+Qh) > Mi {—Qh) in 
Fig. 121 indicates that pinning is stronger when H is closer to the tracks than in the crystal- 
lographically equivalent configuration in the opposite side. This asymmetry demonstrates 
that at the angle +Qh vortices form staircases, with segments trapped into the tracks. Sec- 
ond, for Qh < ©max we again observe asymmetry. Mi (©h) crosses Qh = with positive 
slope, indicating that pinning decreases as H is tilted away from the tracks. Thus, we can 
conclude that staircases extend at least beyond the c axis into the Qh < region. 

Let's now analyze the topology of the vortex staircases. The angle ©^ between the kinks 
and the c axis (see Fig.Ej) can be calculated by minimization of the free energy |^. If Lp is the 
length of a pinned segment, and the length of the kink, the line energy is £" oc LpCp (©d) + 
Lk^f{Qk), where e/(©fc) ~ (©fc) [Iuk + O.S] and epi^Qc) ~ EoE^Qd) [hin + at] are the 
line energy for free and pinned vortices respectively and at < 0.5 parameterizes the core 
pinning energy due to the tracks (smaller at implies stronger pinning). Minimizing E with 
respect to ©^ we obtain the two kink orientations, ©^ for Qh < ©d and ©^ for Qh > ©d- 




FIG. 6: Sketch showing the possible vortex configurations when the field is tilted with respect to 
the columnar defects. 

Since the tracks are tilted, |©^| and |©^[ are different. However, those angles are in- 
dependent of Qh- As |©h — ©dI increases, ©^ remain constant while Lp decreases and 
the number of kinks increases, consequently the pinning energy lowers. This accounts for 
a decreasing Mj as H is tilted away from the tracks. In particular, for Qh = ©^ vortices 
become straight {Lp = 0), thus ip^ = |©^ — ©d| are the trapping angles in both directions. 
In general ©^ must be obtained numerically, but for etan©^ <^ 1 and etan©^) ^ 1 we 
obtain 
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e \ 111 K + 0.5 

Eq.® adequately describes the main features of the asymmetric regionjs^ |6^ in Fig. 
and for = it coincides with the usual estimates^, \^ of ipr- There is, however, 
an important missing ingredient in the standard description presented above, namely the 
existence of twins and Cu-0 layers, which are additional sources of correlated pinning. This 
raises the possibility that vortices may simultaneously adjust to more than one of them, 
forming different types of staircases (see Fig. 

Pinning by twin boundaries is visible in Figure |21 as an additional peak centered at the c 
axis for H = 20 kOe. A zoom in of that peak is shown in Figure [7| The width of this peak, 
~ 10°, is in the typical range of reported trapping angles for twins 0, 0, H, E^- The fact 
that the peak is mounted over an inclined background implies that vortices are also trapped 
by the tracks. Thus, vortices in this angular range contain segments both in the tracks and 
in the twins. These two types of segments are enough to build up the staircases for > 0, 
but for Qh < a third group of inclined kinks with 0^ < must exist in order to have 
vortices parallel to H (see sketches in Fig. |7j). 




-15 -10 -5 5 10 15 



FIG. 7: Zoom in of the data shown in Fig. |21 for H = 20kOe around the c axis. The sketches 
represent vortex staircases for @h > and Qh < 0. 

Another fact to be considered is that the asymmetry in Mj {Qh) disappears as Qh ap- 
proaches the ab planes. This is illustrated in Figure|Hl where Mj data for — \Qh\ was reflected 
along the c axis and superimposed to the results for + \Qh\- There is a well defined angle 
Qsym beyond which Mi (Qh) recovers the symmetry with respect to the c axis (see inset in 
Figure |HI) . 

One possible interpretation is that for 9^ > Qsym staircases disappear, i.e. Qsym = , 
and we are determining the trapping angle (f^ = Qsym — Qd- However, this is inconsistent 
with our experimental results. Indeed, ip^ should decrease with T, and this decrease should 
be particularly strong above the depinning temperature T^p ~ 30 — 40 K due to the reduction 
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FIG. 8: Angular dependence of the irreversible magnetization Mi for three different fields at T = 60 
K. Open symbols: data for @h > 0. Solid symbols: data for Qh < 0, reflected with respect to the 
c axis. Some of the curves were shifted vertically for clarity. The arrows indicate the angle Qsym 
beyond which the behavior is symmetric with respect to the c axis. The procedure of reflection of 
the data is shown in the inset. 

of the pinning energy by entropic smearing effects (section I.C). This expectation is in sharp 
contrast with the observed increase of Qsym with temperature, which is shown in Figure IHl 
for if = 2 T. Thus, the interpretation of Qsym as a measure of the trapping angle is ruled 
out. Moreover, if in a certain angular range vortices were not forming staircases, pinning 
could be described by a scalar disorder strength, then at high fields Mj (Qh) should follow 
the anisotropy scaling lawful'] Mj (H, Q^) = {e (Qh) H). Consistently, we do not observe 
such scaling in any angular range [6^. 

Our alternative interpretation is that, as H approaches the ah planes, the kinks become 
trapped by the intrinsic pinning. This idea has been used by Hardy et al. 43] to explain that 
the Jc at low T in the very anisotropic Bi and Tl compounds with tracks at 0d = 45° was 
the same for H either parallel or normal to them. Our situation is somewhat different, as 
we are comparing two configurations both having kinks. 

We first note that, according to Eq. Q^ cannot be exactly 90° for finite £, thus the 
intrinsic pinning must be incorporated into the model by assigning a lower energy to kinks 
exactly parallel to the ah planes. Vortices may now form structures consisting of segments 
trapped in the columns connected by segments trapped in the ah planes, or alternatively 
an inclined kink may transform into a staircase of smaller kinks connecting segments in the 
planes (see Fig. IH)). We should now compare the energy of the new configurations with that 
containing kinks at angles G^. This is equivalent to figure out whether the kinks at lay 
within the trapping regime for the planes or not. The problem with this analysis is that, as 
G^ are independent of Qh-, one of the two possibilities (either inclined or trapped kinks), 
will be the most favorable for all Qh- Thus, this picture alone cannot explain the crossover 
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FIG. 9: Temperature dependence of @sym defined using the criterium showed in Fig. |H| The sohd 
hne is a guide to the eye. 

from an asymmetric to a symmetric regime in Mj (Qh)- 

The key concept to be considered in this scenario is the dispersion in the pinning energy 
of the tracks |3|5i^]- The angles 9^ depend on the pinning strength of the adjacent tracks {at 
in Eq. ®), thus a dispersion in at imphes a dispersion in 6^. As Qh increases, it becomes 
larger than the smaller 0^'s (that connect the weaker defects) and the corresponding kinks 
disappear. The vortices involved, however, do not become straight, but remain trapped 
by stronger pins connected by longer kinks with larger 0^. This process goes on as Qh 
grows: the weaker tracks progressively become ineffective as the "local" is exceeded, and 
the distribution of 6^ shifts towards the ab planes. When a particular kink falls within 
the trapping angle of the planes, a switch to the pinned-kink structure occurs. In this new 
picture, the gradual crossover to the symmetric regime as \Qh\ increases takes place when 
most of the remaining kinks are pinned by the planes. 

If kinks become locked, the total length of a vortex that is trapped inside columnar 
defects is the total length of a track, independent of G^, and the total length of the kinks 
is oc |tan9/f — tanOol- As |9i^| grows, the relative difference between the line energy in 
both orientations decreases, an effect that is reinforced by the small line energy of the kinks 
in the ab planes. If kinks are not locked but rather form staircases, taking into account that 
the trapping angle for the ab planes is small0, the same argument still applies to a good 
approximation. The temperature dependence of Qsym is now easily explained by a faster 
decrease of the pinning of the ab planes with T as compared to the columnar defects. 

Additional evidence in support of our description comes from transport measurements in 
the dc-flux transformer configuration where, in the liquid phase in twinned YBCO crystals 
(without CD) vortices remain correlated along the c axis for all field orientations^^, in 
contrast to the observed behavior in untwinned crystals. This suggest that, for all angles, 
vortices are composed of segments in the twins and in the ab planes. 
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C. Influence of the Material Anisotropy and Sample Geometry 

In the previous section we showed that the presence of columnar defects in HTSC increases 
the critical current due to the strong pinning and the reduced vortex wandering when flux 
lines are trapped into the pinning potential. We have also shown that at high fields, Jc{Oh) 
exhibits a peak when the applied field is aligned with the tracks orientation. As field 
decreases, two distinct phenomena progressively and simultaneously appear. On one hand, 
we observe a plateau which reflects the existence of a lock-in phase. On the other hand, 
we noticed that the peak shifts towards the c axis. In this section we present a throughout 
study of the origin of this displacement and demonstrate that this shift is a consequence of 
the misalignment between the external and internal field owing to the competition between 
anisotropy and geometry effects. 



1. Relationship between Internal and Applied Fields 

In thermodynamic equilibrium, the internal field B is determined by minimization of the 
free energy|| G(B) = F(B) - £ + ^^"f^^ , where H = B - 47r(l - z>)M. The components 
of the demagnetization tensor at the sample principal axes are {ux, UyjUz), with i^x + i^y + 
Vz = 1. We adopt the notation that z coincides with the crystallographic c axis, and that 
the X axis is perpendicular to both c and H. By standard minimization of G(B) with 
respect to By and B^, and using the free energy F(B) for the intermediate fields regime 
{Hci H Hc2) we obtain, 

. n ^ fi^y^^z^^) sin(29B) In /i + 1 

where fivy, Vz, ^) = (1 ~ ^z) — (1 — ^y)^"^^ is the direction of the internal field, and the 
reduced field h = H/H,2{Qb,T). 

The result ([7j) only assumes uniaxial anisotropy and the coincidence of one principal axis 
with the c axis, and it shows that under those very general conditions the misalignments 
due to both mass anisotropy and sample geometry have the same field and temperature 
dependence. The function /(z>, e), which contains the combined effects of geometry and 
anisotropy, is the key ingredient of the low field behavior, as its sign determines whether b 
leads or lag behind Qh- 

To be more specific, let's consider the typical platelike shape of all the studied single 
crystals, with thickness 6 along the c axis much smaller than the lateral dimensions and 
Ly. To a first approximation = 6/Lx and I'y = S/Ly, thus z/^;, z/y, (1 — z/^) -C 1. If the 
material is strongly anisotropic and the crystal is not too thin, then (1 — z/^) > (1 — i'y)s'^, 
thus / > and 6^ > We will call this the "anisotropy-dominated" situation. In 

contrast, for thin enough samples of a not too anisotropic material (1 — z/^) < [1 — Vy)e^ , 
so Qb < ©H- This is what we will call the "geometry-dominated" case. The extreme 
limit of this case, with an infinite slab (z^x = ^'i/ = 0) and ignoring the anisotropy, has been 



discussed by Klein et a/.|41|. It is also worth to note that for an infinite cylinder with axis 



perpendicular to H, where the geometry effects are expected to cancel out, z/^ = and 
1 

l2' 



= = ^, thus / oc (1 — and Eq. ((7j) reduces to the well known expression for the 



bulk 3, 0; 

The Eq. ((7j) allows us to determine which should be the vortex direction © b for a given 
angle Qh of the controlled variable H. In other words, we might use the columnar defects 
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as detectors of the internal field orientation taking profit from the fact that Jc maximizes 
when Bb = Qd- Thus, if we know Qd and Qmax-, we are able to determine the misalignment 
— = 6d — Qmax- Although the sign of such misalignment is solely determined by 
the sign of /, its magnitude also depends on additional factors such as sin(20£)) and h?. 
Besides that, the misalignment is strongly temperature and field dependent. It is easy to 
see from Eq. ((7j) that Qb ^ for large enough h. 

At this point it is important to note that, although the misalignment between B and H 
is a low field effect, Eq. ((7j) can only be used for H ^ Hd. It turns out that all our data are 
well described by eq. However, the very dilute vortex limit is conceptually interesting, 
and a detailed discussion about it can be found in Ref. 

Table I summarizes the information about geometrical dimensions, mass anisotropy, dose- 
equivalent matching field 5$ and angle 6^ of the CD with respect to the c axis, for all the 
crystals that we will refer in this section. 

TABLE I: Irradiation and shape specifications for all the crystals studied in this section. Crystal B 
was grown at the T.J. Watson Research Center of IBM and irradiated at the Holifield accelerator, 
Oak Ridge (USA) with 580 MeV Sn^°+. 

Crystal material Bq,{kOe) Bd 5{^im) Ly(nm) Lxilim) Uy{xlO~^) ^^(xlO"'^) /(£-, e){ X 10~^) 



A 


YBa2Cus07 


7 


30 


32° 


8.5 


210 


630 


40 


13.5 


+34 


B 


YBa2Cu307 


7 


30 


30° 


20.9 


715 


2150 


29 


9.7 


+19 


C 


YBa2Cu30r 


7 


57 


30° 


11.5 


1050 


1050 


11 


11 


+1.8 


D 


YBa2Cu307 


7 


22 


57° 


4.3 


381 


762 


11.3 


5.6 


-3.2 


E 


NbSe2 


3 


0.5 


27° 


7.7 


765 


640 


10.1 


12 


-67 


F 


NbSe2 


3 


0.5 


27° 


7.7 


585 


640 


13.2 


12 


-63.5 


G 


NbSe2 


3 


0.5 


27° 


7.7 


419 


640 


18.4 


12 


-58 



2. Anisotropy- dominated case 

In Fig. 121 we have already shown the data corresponding to sample A where we noticed 
that the maximum in moves towards the c axis as H decreases. As shown in Fig. 
this effect persist at T = 70 K. A similar behavior is observed in crystal B (see Fig. irUjl . 
although the shift turns out to be smaller than in A at T = 70 K. These two crystals have 
the same anisotropy and irradiation conditions, but different shapes. Thus, at the same T 
and H aU factors in Eq. ^ are identical, except for /(i>, e). As seen in Table I, the difference 
in demagnetizing factors results in a smaller /(z/, e) for sample B than for A. Hence, the 
misalignment in sample B is expected to be smaller, as indeed observed. 

3. Compensated case 

A striking result predicted by Eq. ((7j) is that the competing effects (anisotropy and 
geometry) could be exactly compensated if one were able to tune the demagnetizing factors 
and the anisotropy in order to get fiy^e) = 0, a condition that is satisfied for 1 — z/^ = 
(1 — Uy) e^. For the YBCO single crystals, e ~ 1/7, and this requires extremely thin samples 
with a big area. Table I shows that crystals C and D almost satisfy this compensating 
condition, as the absolute values of their fiy^e) are, respectively, a factor of ~ 20 and ~ 10 
smaller than in A. The right panel of Fig. ITUl shows the angular dependence of M, for these 
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FIG. 10: Left panel: angular dependence of the irreversible magnetization for several fields at 
T = 70 K for the samples A and B (YBCO). The arrows indicate the angular position, Qmax^ 
of the maximum in Mi^Qn) for the lowest fields. Right panel: irreversible magnetization as a 
function of the relative angle Qh — ©d for several fields at T = 60 K in samples C and D (YBCO). 
For clarity, some curves have been translated vertically. 

two crystals at T = 60 K. Since the CD orientations ©£> in samples C and D are different, 
in order to compare them we set the abscise as the relative angle Qh — ©d- In this figure 
we clearly observe that the peaks remain locked at the tracks direction even for the lowest 
fields, in complete agreement with the expectation. The same behavior was observed for 
other temperatures. 



4- Geometry- dominated case 

So far, the samples studied were YBCO crystals with the same anisotropy but different 
geometries. On these samples we observed that the peak either shifts in the direction of the 
c axis or almost does not deviate from the CD direction. As pointed out previously, this 
behavior arises from the strong anisotropy effect in this material. In order to change the 
sign of the deviation {i.e., a shift toward the ab plane), we need to reduce the anisotropy 
effects. (Tabic I shows that crystal D has / < 0, thiis strictly speaking it is in the geometry- 
dominated case, but the shift is too small to be detected). 

To that end we decided to measure NbSe2 single crystals, which have £ ~ 1/3, making 
the anisotropy effect about 5 times smaller than in YBCO. Beside this, very large and thin 
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NbSe2 crystals can be readily found, and they can be tailored at will to obtain the desired 
shape. Thus, we irradiated a rectangular crystal (labelled as sample E), such that it is in 
the extreme geometry-dominated case, with f{u,e) <^ 0. Fig. [TlT a) shows Mj(9j^) for this 
sample at T = 4.4 K for several H. At high fields we observe a peak at the tracks' direction 
0£) = 27°. As field decreases the peak becomes broader and, in contrast to the YBCO 
behavior, it progressively moves away from toward the ah plane, in agreement with a 
negative fiy^e). 




■ 15 30 45 60 75 90 
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FIG. 11: Angular dependence of the irreversible magnetization for sample E (NbSe2) at (a) T = 4.4 
K for several fields (in units of kOe), (b) H = 0.15 kOe for several temperatures. 

The conclusive evidence that the function /(z/, e) dominates the behavior of the misalign- 
ment 0_B — Qh comes from samples F and G, which are pieces of crystal E. These samples 
were obtained by cutting the sample E along a line parallel to its shortest side, in such a 
way that the demagnetizing factor Vx remains unaltered, but Vy increases. In this way, the 
absolute value of the function fiy^e) was progressively reduced, i.e., we moved away from 
the "extreme geometry-dominated case" and approached the "compensated case" (see Table 
I). According to Eq. ((Tj), the deviation of the maximum in Mj(B//) for given H and T should 
become progressively smaller for crystals F and G as compared to crystal E. This is in fact 
observed, as demonstrated in Fig. [T21 where Qmax — ©d for crystals E, F and G is plotted 
as a function oi h = H/Hc2(T, Qh)- 

The misalignments for all the YBCO crystals shown in Figs. |21 and are also included 
in Fig. ^1 Thus, this figure summarizes all the samples studied in the present work, at 
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FIG. 12: Deviation in the maximum of the irreversible magnetization with respect to the tracks' 
direction Gd — © max 1 function of h = H/Hc2, for all samples studied in the present work. 

various temperatures and fields. The three possible low field behaviors are clearly visible: 
anisotropy-dominated (upward curvature), geometry-dominated (downward curvature) and 
compensated (almost horizontal curves). It is worth to note that, in all the not-compensated 
cases and for both materials, the misalignment between B and H becomes relevant for fields 
below a certain characteristic field H ~ 0.02Hc2- 



5. Quantitative test of the model 

The above results clearly demonstrate that the qualitative differences in the low field 
behavior are controlled by the factor f{0,e). We now want to verify whether the H and 
T dependence of the shift is well described by the model. According to Eq. ((7j) these two 
variables appear only through the combination h = H/Hc2{Qb,T). Thus, |0d — ©maxi 
should increase not only with decreasing H at fixed T, as already seen in figs. ITUl and fTlT a) . 
but also with decreasing T at fixed H, due to the increase in Hc2{T). This second expectation 
is also verified, as shown in Fig. ITlT b) where Mj(0//) for sample E was plotted at constant 
field H = 0.15 kOe for several temperatures. 
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The equivalence between the variations in T and H is quantitatively verified in Fig. IT^ 
where sin(9max — ©d)/ f{i^, £) is plotted as a function of h for the two sets of data shown in 
Fig. ITlT a) and lTlT b). We observe a good scaling, thus confirming that h is the appropriate 
variable. The upper critical field values Hc2 = 51.8 kG(l — t) were taken from the literature 
[g^, thus the superposition of the two curves involves no free parameters. 
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FIG. 13: Main panel: sin(Omaa: — 6z))//(i>,e) vs h for the two sets of data shown in Fig. [TTT a) 
and lllf b) for the sample E. Solid symbols: H fixed and T swept. Open symbols: T fixed and H 
swept. Inset: the same scaling shown in the main panel for samples E, F and G. The solid line in 
both, main panel and inset, corresponds to Eq. with Qb = &D = 27° and k = 5.6. 

Finally, we analyze the quantitative effect of the factor /(z>, e). This factor is a constant 
for a given sample, so it is the same for all the data in the main panel of Fig. ^1 In contrast, 
in the inset we show the same scaling procedure for the crystals E, F and G at T = 4.4 K, 
so now /(z>, e) is different for each sample, while all the other parameters remain identical. 
We again obtain a good superposition of the data, although the scaling is poorer than in 
the main panel, probably due to the damage produced in the crystal after each cut process. 

The solid line in the main panel of Fig. depicts the expectation of Eq. ((7j), with 
©B = ©D = 27° as experimentally determined from the location of the maximum at high 
fields, and a single fitting parameter k = 5.6. The same curve is shown in the inset. The 
value K = 5.6 is smaller than the accepted value |^ k ~ 9, a similar discrepancy to that 
observed when studying the lock- in effect. 

6. The Role of the Columnar Defects 

It is important to keep in mind that the shift in Qmax at low h is not due to the CD. 
We are only using them as a passive non perturbative tool to measure the vortex direction 
in the bulk of the samples, what is not easy to do by other methods. In fact, the pinning 
of the CD is not contained in Eq. ((7j), which arises from the minimization of a free energy, 
and thus describes a state of thermodynamic equilibrium. It is obvious, on the other hand, 
that the uniaxial pinning of these correlated structures is relevant and should be included 
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in the analysis. This is usually done by adding to the free energy a term Fpin that accounts 
for the correction to the vortex self-energy due to the CD, and then comparing the energy 
of alternative configurations j3. l43| . 

This additional contribution depends on the orientation of the vortices, Fpj„ = Fpin {^b), 
and it is always negative, reflecting the fact that, for 7^ ©d, the CD promote the 
formation of staircase vortices whose self energy is lower than that of a straight vortex at 
the same average orientation. Fpin (©b) decreases as 9^ approaches Qd due to the increase 
of the core trapped fraction, and it minimizes for that orientation, when the vortex cores are 
totally trapped into the tracks j^]- The key point in the context of the present study, however, 
is that the incorporation of Fpin (©b) into the scenario does not modify the previous results, 
as we show below. 

Let's first consider that H is applied at the angle Qh = ©mas such that, in the absence 
of pinning and according to Eq. ((Tj), the vortices would be at the angle Qb = ^d- If we now 
"turn on" Fpin (©b); the only effect will be to deepen the already existing minimum of the 
free energy at this orientation, without changing the angle. 

Let's now suppose that H is applied at an angle Qh slightly smaller or slightly larger 
than Qmax- In the absence of pinning vortices would respectively orient at angles Qb slightly 
smaller or slightly larger than Qb, according to Eq. ((7j). The addition of the term Fpin (©b) 
will now shift the vortices towards B^, that is, a kind of effective angular attractive potential 
towards the CD orientation will develop. In particular, for \Qh — ©maxl < ^l, the influence 
of Fpin (©s) will be so strong that the system will minimize its free energy by orienting the 
vortices exactly along the CD. In section II. B. we have extensively studied this effect, that 
manifests in our measurements as a plateau in Mj(B//) of width 2(pB. Note that the center 
of the plateau coincides with ©max- Thus, although the relation 9^ vs Qh will be modified 
by the CD, the angle ©max, experimentally defined as the maximum in M((0//) or as the 
center of the plateau where necessary, will still be given by Eq. ([Tj). 



III. AC RESPONSE: DYNAMICS REGIMES WITH COLUMNAR DEFECTS 
NEAR THE SOLID-LIQUID TRANSITION. 

In the previous section we have used dc magnetization measurements to identify and 
characterize vortex structures at different orientations deep in the vortex-solid phase, that 
is, well below the irreversibility line. What changes in this picture if we increase T and 
get very close to Hi„-{T)1 A first issue is that both the vortex characteristic energy Eq and 
the effective pinning energy Ep decrease (the last one particularly fast due to the entropic 
factor), so the lock- in angle tends to zero. As Hc2{T) also decreases, the misalignment 
effects described in II. C, which are proportional to Hc2/H (see Eq. (j7))), also diminish. 
Thus, (except perhaps at extremely low fields) we expect a sharp peak in Jc at the angle 
G15. A point that is far less obvious is whether staircases will still form for all orientations, 
what is their structure, and whether the various correlated pinning mechanisms will remain 
coupled. 

Of course Jc{T) oc Ep tends to vanish (for any orientation of H) as we approach Hirr, thus 
magnetization measurements are not sensitive enough to explore this regime. So we will now 
turn our attention to the ac susceptibility, which has also been extensively used to explore 
the vortex dynamics in HTSC with both correlated an uncorrelated disorder. Indeed, this 
is a natural complement of the dc magnetization studies. The broad range of amplitudes 
and frequencies of the ac field easily accessible with this technique allows the exploration 
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of the dynamics of the vortex system not only in the critical state but in a variety of other 
dynamic regimes. 

When comparing both types of studies, a key factor to be taken into account is the char- 
acteristic size of the vortex displacements involved. When measuring magnetization loops, 
the field changes produce vortex displacements much larger than the distance between de- 
fects. The excitations explored by creep measurements, on the other hand, are half loops, 
double kinks and superkinks (see section I.C), which expand until the vortex (or bundle) is 
completely removed from the CD. So all these measurements involve vortex motion through 
the pinning landscape at an inter-valley scale. In contrast, ac susceptibility measurements 
can be carried out at ac field amplitudes so small that pinned vortices only perform trans- 
verse oscillations that are uniform along the field direction and whose amplitude is only a 
fraction of the pinning range of a columnar defect, i.e., intm-valley motion. This particu- 
larly interesting situation, which results in a linear response with a very scarce dissipation, 
allows to sense the curvature of the pinning potential wells [1, [6^1 . 

In this section we will use ac susceptibility as a tool to first investigate the angular 
regimes that exist in the vortex system in the presence of CD in the vicinity of the solid- 
liquid boundary, and then to explore the dynamic behavior of those regimes in a wide range 
of current densities and excitation frequencies. We construct the dynamic phase diagram 
of the system in the temperature vs. ac magnetic field plane and we determine the various 
crossover fields and currents, as well as the parameters that characterize the different pinning 
regimes. To achieve that, we have developed a couple of useful procedures for the analysis 
of the data. It will become apparent that the characterization of the dynamic regimes is 
essential for the correct interpretation of the angular dependence in the ac response. 



A. Experimental technique 



To measure the complex AC susceptibility x = x' + ^x" i crystals were fixed to one coil of 
a compensated secondary pair rigidly built up inside a long primary coil that produces an ac 
field /iae*^*n which is parallel to the c axis and is very homogeneous in the volume occupied 
by the sample. The amplitude ha can be varied from 5 mOe to 8 Oe and the frequency 
f = ui/2tt from 300 Hz to 100 kHz. A uniform static field Hdc up to 1300 Oe is added. The 
AC coils setup and the sample can be rotated with a precision of ~ 1°, allowing to vary the 
angle G between the c axis and Hdc- (Notation: At the high temperature, near Tbgi where 
these measurements are performed, we can consider B || H. Thus, throughout this section 
we will refer simply to the angle G between the DC field and the c axis, G = G^ = G^. We 
also refer to the angle if = Q — Qd between Hdc and the tracks.) 

Measurements are performed in all the cases by slowly warming up the sample. Generally, 
several curves at various amplitudes ha or at different angles G were recorded during each 
temperature sweep. The scan in angle was always performed in the same direction to avoid 
backlash problems. All the curves of x' and x" at each frequency are normalized by the 
same factor, corresponding to a total step Ax' = 1 with H^c = 0. We carefully looked 
for and confirmed the absence of thermal decoupling between sample and thermometer, as 
well as measurable heating effects due to the excitation current in the primary coil at high 
amplitudes. 

In each sample, we define the critical temperature Tc as the onset of the AC transition 
at Hdc = in the linear AC regime (very low ha), and the zero field transition width ATc as 
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the range with nonzero dissipation in the same condition. The low irradiation doses apphed 
in all the samples reported in this section, produce no measurable effect either in the critical 
temperature or in the transition width. 

B. Angular dependence of the ac response. 

In this section we will identify several angular ranges where pinning is dominated by differ- 
ent mechanisms. With this scope, we will compare the overall angular response of irradiated 
and virgin samples. We present measurements performed on two twinned Y Ba2Cu^0i sin- 
gle crystals from the same batch, both having the same = 91.8 K and ATc ~ 0.6 K. One 
of them (sample A) was irradiated at room temperature with 291 MeV Av?'^^io\is with a 
dose-equivalent matching field ~ 700 Oe at an angle 0_d = 30° from the c axis, and the 
other one (sample B) was used as a reference. 

Figures ^^a) and ll4r b) show some of the experimental x curves recorded as a function 
of temperature with Hdc = 360 Oe at different angles 6 in both the virgin and irradiated 
crystals. The zero field transitions are also included as a reference. Comparison of both 
figures shows that when the direction of Hdc is close to the ah planes (curves for O = 92° 
in Fig. IT^ a) and B = 91° in Fig. IT^ b). shown with solid squares) the response of both 
samples is very similar, indicating that the pinning properties in that configuration are not 
significantly altered by the introduction of CD. In contrast, at all other field orientations the 
irradiation produces a clearly visible upward shift of the X'iT) curves, due to the additional 
pinning introduced by the CD. Close to the defects (6 = 28°, shown with solid up triangles 
in the Fig. Efb)) pinning is drastically increased. 

In the virgin sample, curves at intermediate angles (e.g. O = 24.5°, shown with solid 
up triangles in Fig. Efa)) present a complicated structure. This characteristic structure 
extends to higher B as ha decreases, occurring up to around 75° in the linear regime, as can 
be seen in the inset of Fig. IT^ b). A similar behavior has been reported in recent works 
(67| . Generally, the sharp onset is related with a first order transition and the previous 
decrease in dx IdT is associated with a peak effect jH^, although this last question is still 
open. Moreover, it has become clear by now that measurements in nonirradiated crystals 
at intermediate angles must be done taking into account the thermomagnetic history j67|, so 
data at intermediate angles in Fig. IT^ a) must be taken with caution. In any case, in the 
irradiated sample both the structure and memory effects disappear completely. 

Data as those shown in Fig. EJcan be used to build up curves of x' + 1 as a function of 
B in both samples at fixed T. Figure ITSl shows such curves at three different temperatures. 
Several features are apparent here. First, in the irradiated sample the original symmetry is 
broken (x'(0) 7^ ) by the presence of the tilted defects. This symmetry breaking 

occurs over a very large angular range; the symmetry is only recovered for |B| > Qsym ~ 75°, 
signaled with arrows in the figure. Second, in the asymmetric range |B| < Qgym the screening 
of the ac field is, at any temperature, larger than the one corresponding to the virgin sample. 
Third, beyond |Bsym| the behavior of both samples is very similar. 

The first obvious conclusion that arises from figure El is that pinning in the irradiated 
crystal is dominated by CD over most field orientations. It is apparent that the CD have a 
directional effect {i.e. they act as correlated pinning centers), as they modify the original 
symmetry. This is true even when Hdc is almost perpendicular to the tracks, between —75° 
and —60°. Note also that, even in that angular region, the pinning due to the CD is much 
more efficient than that present in the non-irradiated sample, as indicated by the much 
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larger screening (smaller x')- In the quadrant of positive angles (the right side of the 
figures), where the track direction is included, the vortex pinning in the irradiated sample 
abruptly decreases around the Qsym angle, whereas in the left side, the screening falls down 
at 9 ~ —60° and, not until 15° beyond this angle, at 6 ~ —Qsym, the behavior of the 
non- irradiated sample is recovered. In section III.E. we will show that the angle in which 
the screening falls down (for this case about —60°) is ha dependent, and is a consequence of 
a change in the dynamical regime. 

Figure El also shows that the response is not symmetric with respect to the defects 
direction. This fact can be easily explained since, as we argued in section II, the anisotropic 
character of the material and the presence of natural correlated defects (twin boundaries 
and ab planes) should have an important role in the structure of the vortices. 

The main point to emphasize from these measurements is that the angular region in 
which the CD act as correlated pinning centers is very large. It can also be noticed that, in 
the range of temperature of Fig. UHl Qsym ~ 75° is nearly constant. This behavior is not 
consistent with the existence of a trapping angle ipT beyond which the correlated nature of 
the pinning should disappear, because at these high temperatures a such angle is expected 
to be very narrow and to decrease fast with T. 

We now turn to the situation for field directions close to the ab planes (|6| > 75°). The 
very similar behavior of both samples in this angular region suggests that the main pinning 
source is the same. A related feature is that, in the non-irradiated sample, a qualitative 
change in the linear response is observed very near B ~ 75° (inset of Fig. Efa)). Beyond 
this angle, no structure is observed at any AC field. The above observations could be 
explained by assuming that Qsym indicates the angle beyond which, for both the irradiated 
and virgin samples, the ab planes become the prevailing pinning centers. However, due to 
the strong angular dependence of the AC response in the vicinity of the ab planes, it is 
clear that a similar study with a better resolution should be necessary to confirm the last 
conclusion. 

It is useful at this point to compare and contrast these results with those obtained by DC 
magnetization, described in section II. A clear coincidence is that both sets of results show 
that correlated pinning dominates for all field orientations. The various angular regimes 
are related to different vortex structures arising from the combined effects of CD, twins 
and ab planes, not to crossovers to uncorrelated-pinning-dominated regimes. It is also true 
that, both near the Bose-glass transition and deep into the solid phase, over most field 
orientations vortices form staircases with segments pinned into the CD. The comparison of 
the Qsym obtained by both techniques is particularly interesting. The similarity of the values 
(see figs. IHl 1^ and [T3|) strongly suggests that the underlying physics is similar. On the other 
hand, we should look carefully at the meaning of Qsym in each case. In the dc magnetization 
studies, our interpretation is that for B > Qsym most of the kinks (which connect the vortex 
core segments pinned by the CD) are locked into the ab planes. This implies that intrinsic 
pinning must be important in that angular range. The ac data on Fig. additionally 
suggest that at that high temperature intrinsic pinning may be the dominant pinning source 
in that range, as the ac response of the crystals with and without CD is the same. As this is 
not so clear at lower temperatures, more detailed studies of the interplay of CD and intrinsic 
pinning near Qsym as a function of temperature should be performed. Finally, it is worth 
to mention that, as expected, no evidence of lock-in phase or misalignment between B and 
Hdc is observed at the high temperatures of the ac studies. 

Once the overall angular dependence of the ac response has been determined, we can 
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use the ac susceptibility measurements to explore in detail the various dynamic regimes 
that occur at a given field orientation as a function of T and if^c, so a more microscopic 
information of the vortex-defects interactions can be obtained. Basically, this is achieved by 
changing the amplitude of the vortex oscillations, which are determined by the amplitude 
ha- By performing these studies at a few selected orientations we can identify differences 
among the various angular regions, as will be shown in the next two subsections. 



C. Dynamic regimes of vortices pinned by aligned columnar defects. 

In this section we concentrate in the ac dynamic response when vortices are aligned to the 
columnar defects {i.e. ip = 0). By varying the amplitude and frequency of the ac field, we 
are able to cover a wide range of current densities. At the lowest amplitudes we observe the 
linear response with very low dissipation, characteristic of the oscillation of pinned vortices 
inside the tracks (Campbell regime) IQ, ||66| , while at the high J limit a critical state develops. 
In between, a large non-linear transition regime is observed. We explore both the intra- and 
inter- valley dynamics in the same experiment, and we construct the dynamic phase diagram 
of the system in the T vs. ha plane. We will present results for two irradiated samples. 
One of them is the sample A used in the previous section. The other one (sample C), was 
irradiated with 280 MeV Sb^^+ at an angle of 15° from the c axis to a dose equivalent to 

~ 3000 Oe. 

1. Linear response 

A linear regime is characterized by a x independent of ha. By analyzing both components 
X' and X" "we are able to establish the nature of the regime. The linear ac susceptibility 
X due to vortex motion in a superconductor is determined jlsl. IgqI ItoI ItH by the complex 
penetration depth Xac(T, B,u!) = \r + iXi. The ratio e = Xi/\r, identifies several linear 
regimes. An ohmic response {i.e. a real resistivity) such as flux flow or thermally assisted flux 
flow (TAFF) corresponds to e = 1. A real Xac (^ = 0) represents a non-dissipative response 
with X" = 0, which is almost the case in the Campbell regime where pinned vortices perform 
intra-valley oscillations, although there is always some residual dissipation, thus e -C 1 but 
flnite. 

Our experimental situation is well approximated by a thin disk in a transverse ac magnetic 
fleld where, according to Brandt |48l |. 



Here ip = R6 /2n\\^, where R and 6 are the disk radius and thickness respectively, c„ and 
A„ are tabulated real constants, and the sum arises from the discretization involved in the 
numerical procedure. Expression |S1 allows us to calculate x'{^R:^) cind x"{^R:^)- We can 
numerically invert this function to compute Xr{T) and e{T) from our x{T) data in the 
linear regime [t^. Figure ITBl shows s{T) so obtained in sample A for Hdc = 360 Oe parallel 
to the tracks. We observe that e{T) remains small and approximately constant, e ~ 0.1, 



N 




(8) 



n=l 
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for T < 91. 3-?^. This behavior is characteristic of trapped vortices oscillating inside their 
pinning wells. At T ~ 91.3 K an abrupt depinning occurs, as indicated by the sudden 
increase in e which in a narrow temperature range of ~ 0.2 K grows close io e = 1 that 
corresponds to the ohmic (flux flow) regime. 

It can be seen from Eq. (jSJ that, when e <^ 1 (low dissipation limit), x is solely 
determined by the real part of the penetration depth, Ar. On the other hand, \%{T, B) = 
A|(T) + ^n(T, B), where and Ac are the London and Campbell penetration depths 



respectively [69|, |70[, being 

where is the Labusch parameter that measures the restoring force of the pinning well. 

Thus, measurements of x{T) in the Campbell regime at various Hdc allows us to obtain 
XciT,B ~ Hdc) and aiiT.B ~ Hdc)- Figure IT7I shows X\ vs. Hdc ~ -B at various tem- 
peratures in both samples. In sample C (figure [T7|^a)), B^ ~ 360 Oe without increasing 
the vortex density. It can be seen from figure ITTT a) that, in sample C, A^ oc Hdc ~ -B, 
so Eq. (jn)) implies a field- independent ai. This important result strongly suggests that 
vortices are individually pinned. The last interpretation is well supported by results in the 
sample A f figure ITTT b)): at low dc fields \\ grows linearly with Hdc, but above Hdc ~ 400 
Oe ~ \Bff) there is a faster increase of A|j. with Hdc, indicating that decreases with B. 
We interpret this feature as the crossover to a collective Campbell response. The values 
oi ul obtained in both samples are very similar, demonstrating that the field-independent 
ul indeed characterize the elementary vortex- tracks interaction in a single vortex pinning 
regime. 

We also detect a small dissipative component, frequency independent within our resolu- 
tion, and which is likely to originate in thermally induced jumps between two metastable 
states of similar energy (two level systems) that characterizes a glassy phase^^. We will 
not extend on this topic here, but refer to Ref. for more details. 



2. The onset of nonlinearity 

As ha increases, the response becomes nonlinear. In crystals with aligned CD the crossover 
from linear to nonlinear response is very smooth, so the exact boundary /^^(T) between both 
regimes is difficult to determine and consequently we adopted an experimental criterion 
explained in ref. H^. This boundary signals the appearance of vortex oscillations that are 
large enough to break the linear approximation in the vortex motion equation. This may 
originate either from intra- valley motion (breakdown of the parabolic approximation for the 
pinning potential at large oscillation amplitudes j^,|6^) or from the increase of nonlinear flux 
creep (inter- valley jumps) as J grows [6 9^. In either case, the decrease of the pinning energy 
as T increases implies that h''^ should decrease with T, as indeed observed. 

As the amplitude of vortex oscillations in the linear limit is proportional |69| to the 
local current density J, which is spatially inhomogeneous, the breakdown of the harmonic 
oscillations occur at different ha in different parts of the sample. For a disk in a transverse 
field, the current distribution in the linear regime can be numerically calculated star ting from 
the experimental values of x' (T) and following the procedure described by Brandt '48] . The 
method allows us to estimate the actual current density J{T,p, ha), where p is the distance 
from the center of the disk. Since J is larger at the disk border [p = R), nonlinearities first 
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occur there, when the local current density reaches the value = J(T, R, /i^(T)). Using 

our experimental estimate of the crossover field h}^{T), we can calculate J\T) as shown in 
the inset of Figure ^1 As ha keeps growing the boundary of the Campbell regime moves 
inward into the disk. 

We can also estimate the transverse displacement of the vortices near the edge of the 
sample (those that perform the largest oscillations), u{p ~ i?). In the linear regime, u is 
related with J in the same position through the simple form (1/c) 00 J{p) ~ aLu{p). The 
vortex displacement at the edge of the sample when linearity breaks down, r;, will be 

n = l-^/{T) (10) 

The radius r; measures the range where the pinning potential Sp (r) can be well approx- 
imated by a parabola. Using values for J\T) and aiiT) calculated from the experimental 
data, we obtain r; ~ 50 A in the sample A (irradiated with Au ions), and r/ ^ 25 A in 
the sample B (irradiated with Sn ions). These values have a striking coincidence with the 
observed radius of columnar defects r/) in each caseji^. It is somehow surprising that even 
at these high temperatures where ^ > r^), and so no drastic change in Sp (r) is expected||J] 
at r = r£), the method still provide us information on the real size of the tracks. 



3. The Critical state 



At the highest amplitudes, vortices perform large excursions outside the pinning wells 
and inter-valley motion will dominate over most of the sample. Above a lower boundary 
h'^(T) the response is well described by a Bean critical state ji^. The large intermediate 
region h^iT) < ha < h'^{T) corresponds to a crossover regime where both the intra- valley 
and inter-valley vortex motion contribute significantly to the ac response. 

To define the lower boundary h'^{T) of the critical state regime, we have developed an 
experimental method based in the analysis of the X'iT, ha) data[32|. We will not explain the 
method in detail here but only point out its basic concepts. The Bean Model assumes that 
there is an uniform persistent density current in any region of the sample where inter-valley 
vortex motion occurs. Due to fiux creep effects, that current density J^j is smaller than Jc 
and depends on the frequency of the measurement. This determines a frequency dependent 
Bean length 

AjT,ha) = (11) 

The screening component x' is only a function of A^ divided by some characteristic sample 
dimension. The dependence of x' on this dimensionless variable is shape- dependent and in 
general difficult to calculate. The key idea of the method is to identify data (at fixed 
frequency) at different T and ha that combine to produce the same Ac according to Eq. 
(jllll . The consistency of the resulting function Ac{ha) is checked at each temperature. A 
critical state is established if, at constant T and u, A^ turns out to be proportional to ha- 

Figure ITHT a) shows the result of such analysis at several temperatures in sample A. In 
all cases we observe that Ac oc ha at high amplitudes, thus proving the existence of a 
Bean regime. At low amplitudes there is a systematic deviation from a straight line and Ac 
becomes larger than expected in the Bean regime, indicating the absence of a fully developed 



29 



critical state. The field h'^{T) that marks the onset of the Bean regime, is shown by arrows in 
Figure ITHT a). Note that in the critical state regime dAc/dha oc 1/J^(T), thus the variation 
of the slope of Ac vs. ha with T (which appears as a vertical shift in the log-log plot of 
Figure lIHfa)) allows us to determine the temperature dependence of Jui(T). 

The great advantage of the above procedure is that it allows us to test the existence of 
a Bean regime, to determine its limits in the ha — T plane and to obtain the temperature 
dependence of J{T) regardless of the sample geometry. The function X'(^c) is unique to each 
sample, although, as expected, we found|32|| that the general shape is similar for different 
YBCO crystals with similar aspect ratio. 

To estimate the absolute values of J we need to estimate the absolute values of Ac, and to 
that end we must rely on some geometrical modeling. We use the result [47^ for a thin disk in 
a transverse field, x = 0-5 when Ac/5 = 0.75. In the right axis of Figure fTHT a) we indicate 
the actual values of Ac obtained in adopting such criterion. The corresponding values of 
the persistent current density (identified as JgokHz (T) to emphasize that it is frequency 
dependent) are shown in Figure El The numerical estimates of Ac (T) and J (T) may be 
affected by a systematic error because our sample is not a disk, but the difference will be an 
overall factor of the order of unity and is irrelevant within the context of our analysis. 

Inspection of Figure El shows that the persistent current vanishes very close to the 
temperature where the Campbell regime disappears and e in the linear limit starts to grow. 
The coincidence of both facts identify the Bose glass transition temperature Tbg- In the 
inset of figure El we can see that the temperature dependence of the persistent current 
density JgokHz in the critical state closely follows that of J' (T) , the limiting current of the 
linear regime. This observation implies that the linearity of the intra-valley oscillations in 
the Campbell regime is lost when the current density at the sample perimeter reaches a 
significant fraction (~ 20%) of the current density that flows in the critical state. 

As mentioned above, the dynamics in this limit is strongly affected by flux creep, thus the 
persistent current density J is much smaller than the critical current density Jc ■ Measure- 
ments at different frequencies allow us to determine the temperature and current density 
dependence of the activation energy, U{J,T) in a wide range of J. We have devised |i32i] a 
procedure to extract this information from the frequency dependence of the x{T) curves in 
the critical state. Figure ITST b) shows the main result from this study. There is a character- 
istic glassy behavior U{J) = J^" (see section I.C.). At large current densities U oc J~^'^, 
the exponent progressively decreasing at low J. A dynamic exponent a = 1.5 is typical of 
collective creep of vortex bundles^, Blatter et al. 0] predict this dependence in the 
regime of large bundles, i.e., for transverse bundle dimension larger than A^. (The decrease 
of a at low J could indicate a crossover to a charge density wave (CDW)-type creep regime). 
We also observe a very strong T dependence jsil] . consistent with creep of vortex bundles at 
high temperature. 



4-. Dynamic diagram 

From the above analysis, it is possible to build up a diagram in the {ha,T) plane to 
identify the various dynamic regimes and analyze the crossovers among them. Figure El 
shows a such diagram, for the sample A at H^c = 360 Oe parallel to the CD, and / = 90 
kHz. The response for T < 91.2 K is characteristic of pinned vortices in a glassy phase. 
Below the crossover field h^^{T), we observe a regime of pinned vortices performing intra- 
valley oscillations inside their pinning sites. The ac susceptibility is well characterized by a 
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Campbell penetration depth and a small dissipative component. Above a second line h'^{T) 
(and again for T < 91.2 K), a Bean critical state forms. The persistent current density 
flowing in this regime is determined by thermally activated processes that involve inter- 
valley collective jumps of vortex bundles. In between these two extreme regimes there is an 
intermediate nonlinear region that occupies a large portion of the {ha,T) plane, as h\{T) 
and h'^{T) differ by more than one order of magnitude as seen in the Figure. 

Within the temperature range T ~ 91.2 — 91.3 K both the Campbell and Bean critical 
regimes disappear. This is indicative of a sudden depinning of vortices at a Bose-glass 
transition that occurs in the range 91.2 K < Tbg < 91.3 K. Above 91.3 K, the behavior 
is indicative of a vortex liquid. A linear response is observed for ha < 25 mOe, but the 
dissipation grows rapidly approaching an ohmic response {e = 1) at T~91.5 K (see Figure 
[TBj) . However, at high ha we still observe a nonlinear response, i.e., we cannot access 
experimentally the unpinned (ohmic) liquid region. This is probably because at our low H^c 
the unpinned liquid regime only occurs in a very narrow temperature range below T^. 

Let's now analyze the glassy phase of the dynamic diagram in terms of existing models for 
vortex motion. According to the simplest scenario [62l], the crossover between Campbell and 
Bean critical regimes should occur at a field h*a{T) where the range of field penetration in 
the later (which is proportional to ha) becomes larger than the range of field penetration in 
the former (which is independent of ha). For a longitudinal geometry (e.g., a slab or cylinder 
in a parallel field) this leads to the simple crossover condition Ac{ha, T) ~ Ar(T) but, in our 
case, a meaningful comparison should involve the actual field penetration ranges rather than 
the penetration depths. The practical criterion that we have adopted is to consider that the 
range of field penetration in the linear and critical regimes is the same when x' = ~V2 in 



We can use Eq. (IT^ to estimate the crossover field /i*(T) from Xr(T) and J(T). We 
emphasize that Xr{T) and J{T) are obtained in totally independent ways; the former from 
the linear regime data (section Ill.C.l) and the latter from the high ha data using the 
geometrical procedure described in section 1I1.C.3. The result is shown in Figure IT^ the 
theoretical expectation for the crossover field h*a{T) coincides very well with the line ha{T) 
that indicates the formation of a Bean critical state. This good agreement supports the 
validity of the basic concept, the critical state develops when the Bean penetration range 
becomes longer than the Campbell penetration range. 

We now discuss the origin and nature of the large intermediate regime. The observation 
of a nonlinear response at current densities well below J^qkhzi the current density flowing 
in the Bean critical state, indicates that these nonlinearities are associated with intra- valley 
oscillations. This can be understood by considering the shape of the pinning potential, that 
at these high temperatures is rather shallow. For the present crystal we found in section 
III.C.2 that nonlinearities develop when the amplitude of the oscillations at the sample 
border reaches a value u{p = R) = ri 50 A. We associate ri with the range where the 
pinning potential Ep (r) can be well approximated by a parabola, which is smaller than the 
total pinning range. This picture explains the smooth nature of the crossover from linear to 
nonlinear response at ha{T). 



each case. This crossover should occur when 
to Eq. (fTTj). when 





(12) 



31 



As was discussed in section I.C., in the temperature and field range of our ac mea- 
surements we are deep into tlie collective pinning regime. This is in agreement with the 
conclusion of section III.C.3. However, the linear response in the Campbell regime is well 
described by a field independent Labusch parameter characteristic of single vortex pinning 
(section III.C.l). This apparent contradiction is easily resolved when the size of the vortex 
displacements involved in each case are taken into account. The linear response corresponds 
to oscillations that are uniform along the field direction, whose amplitude (< 50 A) is much 
smaller than the vortex lattice parameter a ~ 2100 A. In this range of displacements the 
energy variation sensed by a vortex due to the elastic interaction with the neighbor vortices 
is negligible as compared to the variation of the pinning energy, and thus the response is 
solely determined by the interaction of a single vortex with the tracks. We note, however, 
that the transverse localization length is much longer than the defect separation, thus each 
vortex is collectively pinned by many defects. The amplitude of the oscillations refers to 
the vortex center of mass, and the effective pinning potential arises from the contribution 
of more than one track. In the critical state regime at high amplitudes, on the other hand, 
the inter-valley vortex excitations involve transverse displacements of the order of the defect 
separation or larger, and consequently the vortex-vortex interactions play a fundamental 
role in the response. 



D. Dynamic regimes in tilted vortices. 

In this section we use the methods described above for aligned CD, to complete the 
study of the angular dependence presented in section III.B with a detailed analysis of the 
AC dynamic regimes at various orientations of Hdc- We will use the same irradiated crystal 
(sample A, B^i = 30°). The angles selected are representative of the various angular regions: 
Hdc forming an angle very small relative to defects (O ~ 34°, (/j ~ 4° ), Hdc along the 
symmetric direction of tracks relative to the c axis (0 ~ —30°, (/? ~ —60°), and Hdc in the 
angular region where the screening fall down (O ~ —70° , ~ 80° ). For comparison, we 
also include the response with aligned vortices (6 ~ 30°, (/? ~ 0° ) and the analogue response 
of the virgin crystal (sample B) for B ~ 30°. In all cases, the response at very low AC fields 
is linear, and it becomes non-linear above a threshold AC field /i|j(B,T). In the following 
subsection we compare the results obtained in the linear regime for each one of the chosen 
angles, and later we analyze the non-linear response. 

1. Linear regime 

We have used the procedure explained in section III.C.l to obtain the linear real pene- 
tration depth and the Labusch constant in the linear regime. To analyze the B dependence 
of ai, we performed measurements for several values of DC fields along the chosen direc- 
tions in the linear regime, and we calculated the corresponding A|j(x')- We have observed 
qualitative differences in the behavior in the various angular ranges. 

In figure Eni various curves of X%{B) for the chosen values of G are compared at T = 90.5 
K. We can see that the qualitative X%{B) behavior observed when B is aligned with the 
tracks still holds if the field is tilted by a few degrees (6 = 34°), but it changes notably 
for larger tilts, and it is also very different in the non-irradiated sample. In the last two 
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cases, no linear dependence in A|. vs. B at low fields is observed, indicating that a Labusch 
constant independent of field does not exist. It can be seen that the Campbell penetration 
length Ac is much shorter when Hdc || CD than for angles far away (4 or 5 times) and than 
for the non-irradiated sample (10 times). For fields higher than the ratio between the 

Campbell lengths in different orientations decreases. 

Furthermore, we observe that the response at = —30° in the irradiated sample is 
different to the one observed at = 30° in all the range of measured fields; this means that 
up to Hdc ~ 2_B(^ the presence of defects is still important. A noteworthy result is that, 
for low DC fields, the responses at —70° and —30° are very similar. For = —30° near 
Hdc = Bip/2, the curve starts to increase more swiftly, and both responses can be clearly 
distinguished. This fact indicates that even at an angle ip ~ —60° from the tracks direction, 
the occupation of defects play an important role. 

The above experimental results can be qualitatively understood as follows. In section 
III.B. we concluded that vortices are partially accommodated in the defects for all the chosen 
angles. The restitutive constant is higher in the vortex segments that are accommodated 
in the columnar defects. For this reason, while the Campbell regime remains individual, 
the average Labusch constant will increase with the length of pinned segments, i.e. when 
the angle ip (relative to defects) decreases. For large ip the situation is more complex: the 
fraction pinned by the tracks is very small (it can be seen that ai for the non- irradiated 
sample at 30° is just half the value observed at —30° in the irradiated sample). Moreover, 
the pinning force for small displacements is comparable to that induced by neighbor vortices 
and is no longer independent of B. 



2. Nonlinear response 

Measurements in all the available range of ha have been performed. The first notable 
observation is that, in the non linear regime, the angular dependence of the AC response 
is much more pronounced than in the linear regime. For large angles relative to the CD, 
as soon as the linear behavior is lost the x{T) curves get wider in temperature and the 
dissipation peaks notably increase. These facts are resumed in figure |^ where the non 
linear behavior at different angles and their comparison with the non-irradiated sample are 
shown. In figure l^lT a). various X'(^a) curves at the same temperature (90.5 K) for different 
angles are plotted, while in figure l^lT b) all the curves included have the same Xr (and 
therefore they were obtained at different temperatures). 

It can be observed that the AC field at which loss of linearity occurs is lower when the 
pinning due to the CD is less efficient: in the non-irradiated sample (solid down triangles 
in the figures), /i^j is almost one order of magnitude lower than in the case of the irradiated 
sample with field aligned to defects (open squares in the figures). When Hdc || CD, departure 
from linearity is smooth, while for large if angles (solid circles and open up triangles in the 
figures) the departure from linearity is much more abrupt. 

The dependence of x on ha for the highest ha is also very different. Whereas in the case 
of aligned columnar defects this function can be explained in the framework of a critical 
state model, for large (p angles this dependence decreases tending to a new linear regime. 
This behavior is still more notable in the non-irradiated sample. Consistently with the last 
observation the maximum of X"(T) approaches the expected values ji^ for the ohmic regime 
(Xmax ~ 0.4:4 for a disk in transversal geometry). The last remark is displayed in figure 
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There we compare experimental points of x" vs. X' + 1 obtained at different angles at 
ha ~ 6.4 Oe, with those calculated for a disk both in the Bean critical state (continuous line 
in the figure) and in the ohmic regime (dashed line). While in the irradiated sample the 
response at G = 30° is similar to the expectation for a critical state, for angles far away from 
the defects it tends to the ohmic behavior. Notice that, as well as Xmax increases, X'(Xmax) 
also tends to the expected value in a linear ohmic regime. 

As another test to close this picture, we tried to check whether a critical state regime 
is established or not, applying the procedure presented in section III.C.3 to data obtained 
at 6 = —70°. A satisfactory solution was not attained, proving that those data are not 
consistent with the critical state model. 

A possible explanation for this behavior is the following: When ip is large, the effective 
critical current density J^^ is much lower than in the case in which Hdc is aligned with the 
tracks. In fact, it is so much reduced that a critical state cannot be established: as soon 
as the linearity is lost, and the vortex displacements become mainly determined by creep 
mechanisms across the pinning centers, the AC field penetration length increases very much, 
and becomes comparable to the flux flow skin depth. In this condition, a nonlinear regime 
in which pinning forces, activated mechanisms and viscous losses contribute significantly to 
the vortex motion is established. 

Finally, another interesting consideration arises from figure |^ as it allows us to re- 
analyze the angular dependence of x(0) previously shown in section III.B, figure IT^ The 
angle beyond which there is an abrupt increase in the AC field penetration in the left side 
of figure Uni (6 ~ —60°), corresponds to the angle for which, at the particular AC field value 
ha = 0.2 Oe, the system changes from a nearly linear regime to a non linear one. This 
fact can be easily verified crossing the curves X'(^a) at ha = 200 mOe in figure |22la). By 
performing the same procedure at lower ha we can observe that, if angular measurements 
were made at lower AC fields, a smooth angular dependence would extend to orientations 
closer to the ab planes. On the other hand, if we perform the procedure for ha = 400 mOe, 
the abrupt increase of AC field penetration will appear before 6 = —30°. For even higher 
values of ha, the region drastically affected by the CD will only involve the peak close to the 
defects, as it has been reported previously 3|. From all this considerations, we point out 



that the apparent angular region notably affected by the CD is ha dependent. Thus, results 
of angular measurements performed using only using AC amplitude, without characterizing 
the dynamic regimes involved, should be interpreted with caution. 



IV. CONCLUSION 

High temperature superconductors with correlated disorder provide a fascinating frame- 
work for the exploration of novel ideas on vortex physics. In addition, our understanding 
of this topic will be essential for the optimization of YBCO superconducting wires for tech- 
nological applications, as vortex motion in those materials is also dominated by correlated 
pinning, arising in that case from fabrication-dependent microstructures. 

We have explored in detail the angular dependent vortex dynamics in type II superconduc- 
tors with aligned columnar defects introduced by irradiation with very energetic heavy-ions. 
We have used dc magnetization measurements deep in the vortex solid phase, and ac sus- 
ceptibility near the solid-liquid transition. We have shown that aligned columnar defects 
are an excellent tool to test models for vortex dynamics, particularly if they are tilted with 
respect to the crystallographic axes, so their effects can be easily distinguished from those 
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arising from mass anisotropy, sample geometry, twin boundaries and intrinsic pinning. This 
allows us, for instance, to use the uniaxial pinning of the columnar defects as a probe to 
determine the orientation of the vortices inside a bulk material, which in general is different 
from the orientation of the applied fields. 

In some aspects we have found an excellent agreement with the theoretical expectations 
of the Bose-glass model. The field dependence of the lock-in angle follows remarkably well 
the 1/H prediction over the whole temperature range of our measurements. In turn, the 
temperature dependence of the lock-in angle gives strong support to the concept of an 
effective pinning energy dominated by the entropic smearing effect. On the other hand, 
both our ac and dc results show that columnar defects produce effective pinning over a wide 
angular range, and that correlated pinning dominates the scenario for all field orientations. 
One consequence of this is the existence of a rich variety of vortex staircases, with segments 
locked into different correlated structures. 

The complexity of the picture is even larger when the high temperature dynamics is taken 
into account. When the magnetic field is aligned with the tracks the ac susceptibility at low 
ac excitations exhibits a linear response arising from vortices oscillating inside the tracks, 
which is characterized by a Labusch parameter independent of the dc field. This indicates 
that the response is dominated by the individual vortex-track interaction, even though the 
critical current in the same temperature and field range is determined by collective pinning 
mechanisms. The solution of this paradox is related to the fact that the characteristic size of 
the vortex displacements in each case is very different. Finally, the ac response at different 
angles indicates that the characteristics of the dynamic regimes and their extension in the 
temperature vs ac field plane are strongly dependent on the orientation of the vortices with 
respect to both the columnar defects and the crystalline structure. 
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FIG. 14: a) Temperature dependence of the out-of-phase component x of the ac susceptibihty with 
H(jc oriented at various angles relative to the c axis in a nonirradiated sample. Inset: curves of 
X" as a function of T in the linear regime; a structure is present at intermediate angles, b) Same 
in a sample irradiated with a matching field i?^ = 700 Oe at 6 = 30°. 
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FIG. 15: Angular dependence of x' + 1 in an irradiated (open squares) and virgin sample (circles) 
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FIG. 17: Field dependence of the square of the dimensionless hnear real penetration depth at vari- 
ous temperatures in two irradiated samples with different doses. If B < a linear dependence 
holds. 
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FIG. 18: (a) Test that proves the existence of a critical regime. The Bean length Ag is proportional 
to ha above h^iT) signaled by the arrows, (b) Current density dependence of the activation energy. 
The U values came from the frequency dependence of x' in the critical state. 
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FIG. 20: Field dependence of the square of thejlimensionless linear real penetration depth in an 
irradiated sample at a fixed temperature for different angles O. A curve of a virgin sample is also 
shown. 
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FIG. 21: Linear and nonlinear behavior at different angles as a function of the ac field, (a) 
Comparison of X'(^a) curves at the same temperature, (b) Comparison at the same linear real 
penetration depth. Symbols correspond to both panels. 
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FIG. 22: Experimental curves X"ix') ^ large ha = 6.4 Oe in the irradiated an virgin sample for 
various angles. Lines are the theoretical calculation for a disk in an Ohmic regime (dashed line) 
and in a Bean critical state (solid line). 
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